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Abstract 

Calculations of nuclear spin-lattice relaxation rates are carried out by means of exact diagonal- 
ization on small (2 x 6) antiferromagnetic Heisenberg ladders, using the simplest forms permitted 
by symmetry for the hyperfine couplings for the three nuclear sites in CU2O3 ladders. Several val- 
ues of the rung/chain exchange ratio J±/ J\\ have been considered. Comparisons with experimental 
results, field theoretic calculations, and the Gaussian approximation highlight some open problems. 

PACS numbers: 75.10.Jm, 75.40. Gb, 75.40.Mg, 76.60.-k 
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I. INTRODUCTION 



Spin-ladder systems, in particular the two-leg, S = 1/2, antiferromagnetic variety, have 
been the subject of considerable theoretical and experimental investigation.0 Spin ladders are 
appealing because they are one-dimensional systems and thus can be effectively investigated 
using many powerful theoretical tools, while offering a wider parameter space of "simple," 
and potentially experimentally realizable, Heisenberg Hamiltonians than spin chains. The 
simplest Heisenberg ladder Hamiltonian has the form 

'H = ^||('S'n,l ■ Sn+1,1 + Sn,2 " Sn+1,2) + J±Sn,l " 'S'„^2 (l) 

n 

which offers a dimensionless parameter J±/J\\ which is in principle tunable by chemistry or 
pressure. In addition, compounds containing weakly coupled CU2O3 ladders are appealing 
because of possible connections with cuprate superconductivity. 

The present work was motivated by the nuclear spin-lattice relaxation measurements in 
La6Ca8Cu2404i, an undoped ladder compound, by Imai et a/.i These measurements were 
carried out for all of the nuclear sites on the ladder, namely the copper, the "rung" oxygen, 
and the "ladder" (or "chain") oxygen, over a wide temperature range, from low temperatures 
up to nearly 900 K. Because the principal exchange interactions in cuprates are so large, 
on the order of 1000 K, it is quite challenging to do experimental work at temperatures 
significantly greater than the spin gap (A ^ 500 K). 

The experimental results (see Figure 1(c) of Imai et have the following noteworthy 
features. At temperatures below about 425 K, the relaxation rates for all three sites follow 
a common (activated) temperature dependence up to a scale factor. However, on increasing 
T the copper 1/Ti (which we will refer to as l/^"Ti) exhibits a rather sharp departure 
from that of the two oxygen sites (l/*-*(^^Ti and l/'^^'^^Ti for ladder and rung, respectively). 
There seems to be a nearly discontinuous decrease in the derivative of l/^"Ti; moreover, 
above 425 K the if^^Ti data appear nearly linear with an almost vanishing intercept. The 
relaxation rates for the two oxygen sites, in contrast, exhibit no particular features in the 
vicinity of 425 K. 

Several aspects of the wavevector dependence of the low-frequency spin susceptibility can 
be gleaned directly from the data. 

One can express the spin-lattice relaxation rate in terms of the dynamic structure factor 
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for the Cu^^ spins 

^ oc y dqHn{q)S{q,uJn) (2) 

where if„ is the hyperfine form factor associated with nucleus n, uOn is the NMR frequency 
(which we will take to be zero in everything that follows), and S is the structure factor. The 
proportionality constants can be neglected for our present purposes. The spin correlations 
are isotropic, so there is no need to consider the various components, S^^ and so forth, 
individually. The hyperfine interactions are not isotropic, so the orientation of the magnetic 
field in the NMR experiment does affect the results; however, all of the results of present 
interest can be obtained with a single field orientation, which then specifies Hn{q) uniquely. 
The largest hyperfine couplings are between a given nuclear site and the closest spins; at 
that level of approximation, and taking the intra- and inter-chain lattice constants to be of 
unit length, one has 

Hcn = A\ Hon) = 4C2cos2(gx/2), 

Hoi2)=^F^cos^{qy/2) + D^ (3) 

where C, F, and D are the hyperfine couplings identified in Fig. 1(a) of Imai et a/.,i A is 
the on-site hyperfine interaction for copper, and we have elided the orientation dependence 
of the hyperfine interactions (so, for example, A"^ should really be A^ + Ay if the static field 
is along the z axis). 

The essential difference between copper and oxygen sites is that in the latter the hyperfine 
interaction in the vicinity of q = (vr, tt) is much smaller than in the vicinity of q = (0, 0). If, 
at all temperatures of experimental relevance, S{q, 0) had most of its weight in the vicinity 
of qr = (0,0), then then all three relaxation rates would have tracked one another. The 
marked decrease of 1 /'^'^Ti relative to the other two relaxation rates at 425 K indicates that 
this cannot be the case, and in fact suggests that at temperatures below 425 K the ratio of 
the spectral weight near (tt, vr) to that near (0, 0) is roughly constant and of order unity, 
while above 425 K the ratio falls markedly. (The decrease is crucial. If there were an increase 
in 1 /'~^^Ti relative to the oxygen rates with increasing T, one could ascribe that to a turn-on 
of ^((vr, vr), 0) for T > A but ^((Tr, vr), 0) might have been negligible compared to 5'((0, 0), 0) 
at lower temperatures.) 

Why the emphasis on g = (0, 0) and (vr, vr)? In gapped systems such as spin ladders, the 
low-energy spin fluctuations are Raman processes, and at low temperatures one needs to 
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consider only the lowest energy magnons, namely those near q = (vr,7r). Spin fluctuations 
near (0, 0) are associated with two-magnon processes, and those near (vr, tt) with three- 
magnon processes, and on the face of it one would be justifled in neglecting the three- 
magnon processes entirely at low temperatures: see Ref. ^ and references cited therein. 
However, as we have just seen, this appears to be inconsistent with the experimental data 
for La6Ca8Cu2404i, and it is also inconsistent with the quantum Monte Carlo calculations of 
spin-lattice relaxation in a particular Heisenberg ladder {J±/J\\ = 1) by Sandvik, Dagotto, 
and Scalapino,@ at least at temperatures greater than half the magnon gap. 

An extensive theoretical treatment of spin dynamics in gapped one-dimensional Heisen- 
berg models, including spin ladders, has been presented by Damle and Sachdev.i Their 
analysis of S{q, ~ 0) was restricted to q near (0, 0), but they did find the quite interesting 
result that the activation energy for 1/Ti is larger, by a factor of 3/2, than the activation 
energy for the uniform static susceptibility (which is simply the spin gap). An analysis of 
S{q,uj ^ 0) for q near (vr,7r), for systems with J\\ ^ J± has been presented by Ivanov and 
Lee.i Their results are suggestive of a fairly sharp crossover from low- to high-temperature 
regimes at T ^ A, and also indicate that the (vr, tt) contribution to 1/Ti "overshoots" its 
T = oo value and thus decreases as T cxd. 

In the present work, we have applied exact diagonalization to evaluate spin-lattice relax- 
ation rates, following the method of Sokol, Gagliano and Bacci.i We have considered three 
different ladder Hamiltonians, namely J±/J\\ = 0.5, 1.0, and 2.0, and have obtained 1/Ti 
for Cu, 0(1), and 0(2) sites taking the simplest conceivable hyperfine couplings, namely 
A = C = F = 1, with all other interactions neglected. All of the calculations were for rather 
small systems, 2x6, such that exact diagonalization could be carried out in an extremely 
straightforward manner. 

It was noted above that calculations of spin lattice relaxation rates for spin ladders 
have already been carried out by means of large scale quantum Monte Carlo,@ but those 
calculations were limited to the Cu sites. Our goal is somewhat different than that of 
Sandvik, Dagotto, and Scalapino's work. We are not trying to fit the data in detail, rather 
we want to see what can be learned from modest numerical calculations. One reason not to fit 
the data is that to get the gap correct to 10% by exact diagonalization for J±/ J\\ = 0.5 would 
require a system at least 2x12. Another is that we do not treat the spin diffusion contribution 
to the relaxation rates correctly: our calculations effectively introduce an artificial cut-off 
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so that we obtain a finite spin-lattice relaxation rate. Finally, the precise form of the 
spin Hamiltonian for the cuprate ladder compounds is still subject to argument. Although 
the Knight-shift results of Imai et a/.l appear to be consistent with the simple spin-ladder 
Hamiltonian of Eq. (|I|) for J±/J\\ ~ 0.5, it has been suggested by Brehmer et a/.0 that 
instead J±/J\\ ~ 1 and in addition there is a modest amount of plaquette "ring exchange" 
in the Hamiltonian. A quantum-chemical analysis of the exchange interactions in various 
cupratesi provides some support for the latter proposal, since it concludes that J±/ J\\ ~ 
0.9 for Sri4Cu2404i (which is a lightly-self-doped version of the undoped La6Ca8Cu2404i 
compound). 

To be precise, the goals of our calculation are as follows. First, we want to verify that 
S{q,0) has significant weight near q = (7r,7r) as well as near (0,0) and see if there are any 
noticeable trends with varying J±/J\\. Second, we want to explore the crossover from low 
to high temperature behavior in l/Ti. can we see anything like the experimental results, or 
like the theoretical results of Ivanov and Lee? Third, we want to keep our eyes open for any 
unanticipated patterns that might emerge in the numerical results. 

II. METHOD OF CALCULATION AND RESULTS 

The finite-size calculations of spin-lattice relaxation rates are carried out following Sokol, 
Gagliano, and Bacci.i Rather than repeating their discussion of the method let us make a 
few remarks. We take Jy as the unit of energy. 

The first step in the calculation is a complete diagonalization of the Hamiltonian and 
evaluation of matrix elements for certain local spin operators (depending on which nuclear 
site one is interested in). For the 2x6 lattices all of the calculations could be done using the 
simplest possible representations of the states in terms of local values; it was not even 
necessary to use translational invariance to classify states by wave vector. 

The second step is the construction of an auxiliary function which Sokol et al. refer to 
as This is implicitly dependent on T and the hyperfine couplings. We considered 

temperatures ranging from 0.3 to 50. Typically we constructed I{uj) at intervals of 0.02 in 
uj up to at least uj = 0.6. 

Finally, one needs to estimate the zero-frequency derivative of I{uj), because 1/Ti is is 
proportional to T{dl /duj)\^=o. At high temperatures is quite smooth, but at temper- 
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atures comparable to the gap significant structure develops (see Fig. |l|). In order to avoid 
introducing spurious temperature dependences into 1/Ti it is important to use a consistent 
procedure for extracting the derivative from the data. What we did was to fit a zero-intercept 
line through all the data points up to a cutoff cj^nax; weighting all points equally in the fit. 
We did all of the calculations using both cUmax = 0.5 and 0.3. While there are noticeable 
differences in the results using these two cutoffs, as shown in Fig. 0, our conclusions turn out 
the same no matter which is chosen. The use of a much smaller cutoff, which might seem to 
be preferred on the grounds that one is really looking for a zero-frequency derivative, is not 
beneficial. The structure that develops in I{uj) as T is lowered, making it look like a Devil's 
staircase, is a finite-size artifact and must be averaged over, using a suitably large cUmax, to 
obtain results that are representative of the thermodynamic limit. 

We now turn to the results of the calculations for the three nuclear sites and three values 
of Jj_ considered (0.5, 1.0, and 2.0). In every case we take cUmax = 0.5. In Fig. |^ we present 
results on a linear temperature scale, for T < 2. The behavior of the spin-lattice relaxation 
rate at high temperatures is a bit surprising: comparing the plots in Fig. |^(a) through (c) 
it is apparent that while the Cu and 0(1) rates decrease strongly as J± increases, the trend 
for the 0(2) rate is different. This is made more explicit in Fig. ^, where we show 1/Ti for 
all three sites as a function of J_l at T = 50 (effectively infinite temperature). In contrast, 
at low temperatures 1/Ti decreases with increasing J± at all sites, as one would expect since 
the spin gap is an increasing function of J±. 



III. DISCUSSION AND CONCLUSIONS 



It is evident that for Jj_ = 0.5 and 1.0, 1/Ti for all three sites is nearly equal for tem- 
peratures below the spin gap. (Of course we do not claim that this holds to arbitrarily 
low temperatures, just that it seems correct for temperatures as low as we dare to estimate 
1/Ti.) Because of our choice of hyperfine interactions, this suggests that in such cases the 
weight in S{q,0) for q ^ (vr, tt) is approximately three times that for q ^ (0,0). This is 
in quantitative agreement with the results of Sandvik et a/.@ at J± = 1.0. However, the 
story is rather different at J± = 2.0, where the spin-lattice relaxation rates for all three 
sites, including the two oxygen sites, are significantly different even at T = A/2. In the 
strong-coupling limit, then, the simple picture for S{q, 0) in which its weight is concentrated 
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at (0, 0) and (vr, vr) does not work even for temperatures that are a modest fraction of A. 

What can we say about the low-to-high temperature crossover in the spin-lattice relax- 
ation rates? First of all, the sort of behavior seen experimentally, in which 1/Ti for the 
oxygen sites track each other closely while l/^"Ti splits off, appears to be a special feature 
of ~ 1 in the present calculations; it is not at all generic and does not hold for the 
putative experimental value J±_ ~ 0.5. Second, in no case does l/'^^Ti exhibit any sort of 
sharp "break" as seen experimentally; nor does l/'~^^Ti exhibit linear-in-T behavior (with 
zero intercept, or otherwise) in the high temperature regime, even over a restricted tempera- 
ture range (say A to 2A). Finally, in no case does l/'"'^Ti exhibit an "overshoot" during the 
crossover: the spin-lattice relaxation rate associated with all sites monotonically increases 
with T. 

Our calculations thus suggest that there are quite a few open problems in this field. 
Almost none of the prominent experimental facts concerning l/'~^^Ti in La6Ca8Cu2404i are 
reproduced in our finite-size calculations. Furthermore, the work of Ivanov and Leel does 
not seem to have much to say about our results, either. Their calculation is controlled only 
in the J_l ^ 1 regime, so we should only look at the = 0.5 data. Here we have no 
evidence of overshoot in l/^"Ti, and no reason to believe that one can just examine the 
spectral weight near (vr, vr) since 1/*^*^^)Ti "peels off" from l/'-'^^^^Ti in a manner not very 
different from 

At this point we face several alternatives. It is possible that our results are simply 
unreliable, because we are considering systems that are too small (especially for = 0.5) 
and our procedure for estimating dl{uj)/duj is flawed. We cannot rule this out, but we 
strongly suspect that the trends in the results as a function of Jj_ are robust. It is possible 
that the spin Hamiltonian for the ladders in La6Ca8Cu2404i is more complicated than the 
model we have considered. Whether the Hamiltonian of Brehmer et a/.i can reproduce the 
spin-lattice relaxation data requires another calculation. Another possibility that must be 
considered, given the remarkably sharp feature in 1 /^^Ti found in the experimental data, is 
that La6Ca8Cu2404i undergoes, by coincidence, a subtle structural transition at 425 K. This 
could introduce an anomalously strong T-dependence to the hyperfine interactions, though 
why the effect should be so much stronger in Hc^{q) than ifo(i)(9) and i^o(2)(9) is difficult 
to envision. 

Let us now turn to the results of our calculations for spin-lattice relaxation at very 
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high temperatures, shown in Fig. ^. The most natural way to think about these results is in 



/ dqHn{q)S{q, u) is a Gaussian function of u, and then evaluate the frequency cumulants of 
this function by means of short-time expansions of time-dependent correlation functions. At 
T = oo the calculations are especially simple, because the expectation values of correlators 
{Si ■ Sj) vanish for sites i j. For three sites of interest in Heisenberg ladders, the Gaussian 
approximation yields the following exchange dependences of the spin-lattice relaxation rates 
at T = oo: 



can be derived in another way, by considering the strong- J_|_ limit. Then one most naturally 
thinks about the states in terms of singlets and triplets on the rungs. The relevant energy 
scale for the dynamics of the total spin on a rung, which is relevant to l/'-'^^^^Ti, would seem 
to be proportional to Jy (that is, the bandwidth in lowest-order perturbation theory for a 
triplet excitation in a single background^) , and with the hypothesis of a single energy scale 
in f dqHo(2)S{q,uj) one reproduces the Gaussian approximation result. 

We see in Fig. | that 1/Ti for the copper and ladder oxygen sites decreases with increasing 
J_L, qualitatively in agreement with the Gaussian approximation, although the dependence 
on J_|_ is not as strong as that approximation suggests. Furthermore, l/'-'^^^^Ti exhibits an 
increase with J^. The rather poor performance of the Gaussian approximation is somewhat 
disappointing, considering how well it works for estimating spin-lattice relaxation rates in 
square-lattice Heisenberg antiferromagnets.iJi^ It is not too surprising, perhaps, given that 
the dynamic correlations in the S = 1/2 Heisenberg chain are far from Gaussian at T = oo.0 
So, there is yet another open problem in the area of low-energy spin dynamics of Heisenberg 
ladders. 
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FIG. 1: The auxiliary function I{io) for the Cu site with J± = 1.0 and T = 0.5 (squares), 1.0 
(circles), and 5.0 (diamonds). 
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FIG. 2: Estimated values of for J± = 1.0 as a function of temperature, taking Wmax = 0.3 

(+) and 0.5 (x). 
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FIG. 3: Spin lattice relaxation rates as a function of temperature (in units of Jy) for copper 
(circles), ladder oxygen (squares) and rung oxygen (diamonds) sites, for J±/ J\\ = 0.5, 1.0, and 2.0 
in (a) , (b) , and (c) respectively. The upside-down triangle on each graph indicates the value of the 
spin gap A for the corresponding system. 
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FIG. 4: Spin lattice relaxation rates at T = 50 as a function of J± 
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